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ABSTRACT 



A new model for lamellar surfaces formed by anisotropic molecules is proposed. The 
molecules have internal degrees of freedom, associated with their flexible section of length 
N at zero temperature. We obtain a 2D non-standard six vertex model, which is exactly 
soluble and exhibits a finite order transition. The order and the character of the transition 
are determined by the dominant term in the -^-expansion of the interaction energy. The 
dependence of the critical temperatures on N is, instead, determined by the non-leading 
terms in the same expansion. 



1. INTRODUCTION. 



In this work we shall study the statistical properties of a hamiltonian model for a two 
dimensional object formed by flexible molecules arranged in the ordinary 3D space. 

We think that this ideal model not only can mimic real physical surfaces such as layers 
of molecules forming films and interfaces or (semi)bilayers in biological membranes, but 
also can be useful for the interpretation of the layered structures, and related transitions, 
observed in smectic phases of liquid crystals. 

In view of our purposed description, two properties of the elementary components are 
of main interest. The first one is the presence in the molecule of a (long) flexible section, or 
tail, typically an alkylic sequence; the conformational flexibility is associated with internal 
degrees of freedom which represent the rotations around skeleton bonds in the tail. In 
agreement with the formalism of Volkenstein and Flory^ 1 ' 2 ) who developed the conforma- 
tional theory of chain molecules, we shall refer in this paper to trans planar chains and 
to chain defects such as kinks, antikinks and their combinations; we shall also make use of 
Flory's definition of conformational (intramolecular) energies The second structural 
property of interest is the presence in the molecule of a rigid core, or head, earring an asym- 
metric charge distribution, equivalent to a dipole; the relative orientation of such a dipole 
with respect to the molecular longitudinal axis is also meaningful in our analysis. It is well 
known from experiments that systems formed by molecules with these properties exhibit a 
rich variety of mesophases (first classified by Friedel^ 3 -') which, beyond the perfect crystal, 
go toward the isotropic liquid trough a progressive loss of intermolecular order^ 4-10 ). Both 



thermotropic and lyotropic behaviours have been observed. In thermotropic systems, one 
has a pure substance formed by mesogenic molecules; the transitions between mesophases 
occur by varying the temperature of the sample. The order parameters describe the orien- 
tational and the positional arrangement of molecules; they change with temperature and 
are discontinuous at critical temperatures. In lyotropic systems, one has a solution where 
the solute is a mesogenic molecule; the phase transitions are achieved by changing the 
concentration of the solute or by changing the character (ionic strength) of the solvent. 

In this paper, we shall refer to thermotropic compounds. In this field, a wide col- 
lection of experiments (see for instance ref.10), investigated the relations between struc- 
tural properties of molecules and liquid crystalline basic properties; the main interest has 
been focussed on the spatial organization of molecules in the various mesophases, the 
behaviour of critical temperatures in series of homologues and the sequences of occur- 
ring mesophases. In particular, the layered, head to tail, arrangement of molecules was 
observed, since 1975, by Leadbetter et alS 11 ' 12 ' in the smectic A phase. By X-ray diffrac- 
tion on cyanobiphenyls compounds, they confirmed that the layer spacing is greater than 
the molecular length, consistently with the proposed setting of pairs of molecules, forming 
dimers, with heads facing one an other in the center of the layer and tails pointing outward 
in opposite directions. Analogous layered arrangements were observed in smectic phases for 
many other thermotropic compounds; as examples we may quote the nitroderivatives^ 13 ) 
and cyanatophenylbenzoates^ 14 ) with alkylic tails. In these substances many different Sa 
phases are recognized by X-ray diffraction, calorimetric^ 13 -* and dielectric*- 14 ) measures; in 
each phase the head to tail setting occurs with a different extent of molecular overlapping 



and layer thickness. A great variety of Smectic- Smectic transitions may occur; in particu- 
lar, the transition Sai — > Sa2, from the phase with layer thickness equal to the molecular 
length, to the phase with layer thickness twice the molecular length^ 13 ) is observed for 
increasing temperature; the transition Sai to nematic is also observed* 13 ' 14 ). 

Two notes are in order about these facts. First, the persistence of layered structures in 
the different smectic phases suggests the interaction between adjacent layers to be weak, 
compared with the intralayer forces; for this reason we think that it is meningful to study 
the single, smectic-like layer as a two dimensional, independent, neutral and thermotropic 
object in the 3D space. Second, the head to tail setting of adjacent molecules suggests the 
antiferroelectic arrangement of molecular dipoles, parallel to the longitudinal molecular 
axis, to be compatible with the alignment of molecular axis. 

From a theoretical point of view, the formation of mesophases is equivalent to a 
sequence of symmetry breakings. Starting from the completely symmetric 3D phase 
(isotropic liquid), one expects the strongly asymmetric interactions between elementary 
constituents to promote, step by step, the orientational order, the formation of stable and 
weakly interacting objects with lower dimension than that of the environment (such as 2D 
smectic layers, or ID columnar structures), up to the perfect crystal^ 15 ). However, this 
progam is rather difficult to execute. 

The first step has been overcome by Onsager, since 1949^ 16 ' 17 ^, with an exactly soluble 
model for a system of rigid rods with first order isotropic to nematic (I-N) transition. 
There is a wide and rich literature about this transition; one has the phenomenological 
description by Mayer Saupe* 18 ), the Landau theory of orientational order parameter by 



De Gennes ( 19 \ and the order parameters by McMillan's model^ 20 ' 21 ). The presence of 
a flexible section in the molecule, and its effects in the I-N transition, was considered in 
Marcelja's self consistent analysis of chain ordering in liquid crystals^ 23 ) and in Pink's model 
of odd-even effect^ 24 ); the problem was also investigated by Martire^ 25 ) and Luckhurst^ 26 ) 
with molecular field methods. The orientational transition was also studied by Monte 
Carlo simulations; we recall the first calculation by Lebwohl and Lasher ( 27 ) and, more 
recently, the planar versions of the RP n model^ 28 " 30 ). 

For other transitions, between different smectic phases (S-S) or for smectic to nematic 
transition (S-N), the research is not so widely developed; in particular, there are no ex- 
actly solved models and one turns to phenomenological theories of order parameters or to 
numerical simulations. For liquid crystals, the translational order parameter for S-N tran- 
sition was defined by McMillan^ 20 ' 21 ) and compared with the phenomenology of smectic 
phases^ 22 ). The problem about the role played by conformational flexibility in these phases 
and related transitions is avoided: the molecules are ever considered as rigid rods. 

On the other hand, in the field of biological membranes, Nagle^ 31 ^ and, independently, 
Marcelja ( 32 ) developed two models for the transitional behaviour of molecular chains in 
smectic-like phases of lipid mono and bilayers. In these models, the attention is pointed on 
the conformational collapse of the constituents rather than on the effect of such a collapse 
on the surface of the biomembrane. 

Before introducing the model, which is addressed to the study of S — S and S — N 
transitions, we note that the criticality of layered structures is of main interest in a rather 
different context, namely in studying the interface formation between (e.g) two solids or 



film deposition; in this case the geometry of the surface is characterized by the formation 
of ledges and terraces, with the typical roughening transition. Surface structures and 
related phase transitions have been studied with exactly soluble models^ 33 ); we recall the 
ones by Van Beijeren (BC-Solid-On-Solid ( 34 > 35 )) and by Burton (Terrace-Ledge-Kink ( 36 )), 
firstly introduced for the crystal-vacuum interface. The common feature of these models 
is that they can implement a conservation law by a six vertex lattice system, with exactly 
known ( 37 ~ 39 ) solution. Monte Carlo calculations on specific realizations of BCSOS or 
TLK models simulate (e.g) adatoms forming thin films ( 40 > 41 ) or interfacial molecules ( 42 \ 

In this work, we present an exactly soluble model for a single molecular surface, namely 
for a 2D object like a free smectic layer of a thermotropic liquid crystal. Together with 
the general structural properties previously required for the constituent molecules (the 
presence of one head and one tail) , we consider the following distinctive features: 

a) - we are interested in a system whose mtramolecular energies, associated with 
conformational chain defects, are of the same order of magnitude as the mtermolecular 
energies; 

b) - we represent the charge distribution associated with the molecular head by a 
dipole parallel to the longitudinal axis of the molecule; 

c) - we assume an antiferromagnetic structure of the foundamental state; the con- 
stituents are aligned with an head to tail setting. 

There are at least two physical reasons for assuming the property a). First, in many 
molecular systems with smectic transitions, the transition heat is of the same order as 
the melting heat of the crystal formed by the tails only (typically n-alkanes lamellae) ( 4 ) . 



Second, in the same systems, IR and Raman spectra directly show that the alkylic part 
of the molecular layer is liquid like at the critical temperature; in fact, the spectroscopic 
markers of the conformational defects characterizing the liquid phase of n-alkanes are 
clearly observed ( 43 ^. 

As for the requests in b) and c), we stress that they are not independent one each other. 
In fact, constituents with head dipoles orthogonal to the longitudinal axis, should be 
compatible with a layered structure in which molecules have an head to head setting (like 
in a semi-bilayer, in which the alkylic chains point in the same direction). In this case we 
expect to be addressed to the phenomenology of lyotropic systems, in which the transitional 
properties are determined by the dominant interaction between polar heads and solvent; 
alter naltively, for thermotropic compounds, we expect the system to have a transitional 
behaviour dominated by vortex formation in the plane of molecular heads^ 44 ). 

Assuming longitudinal head dipoles instead, the arrangement of molecules, in a temper- 
ature regime with high orientational order, is compatible with the head to tail setting. 
We think that, in this case, there is an effective competition between the dipole-dipole 
and the pure chain-chain contributions to interaction energy; pairs of opposite dipoles will 
have a tendency to reach an antiparallel setting with aligned centers; on the other hand 
tails will prefer to attain the full extended trans planar conformation with aligned centers 
of mass. In presence of a rigid core, this competition leads to an interaction potential 
with (at least) two different minima which, roughly speaking, correspond to the maximum 
coupling of dipoles and chains respectively. Let us consider, in the 3D space, the packing 
of pair of molecules in one of the equilibrium relative positions; we shall obtain the lay- 



ered smectic-like structure, namely the single surface of interest. Furthermore, different 
geometries select different classes of conformational defects, which are compatible with 
alternative packings. We can study the stability, the thickness and the degree of confor- 
mational disorder of such a surface as a function of both the conformational energy and 
the energy difference between the minima of the two-body potential. 

In section 2. we construct the bare model of the surface, obtained neglecting the con- 
formational degrees of freedom. We show that it belongs to the class of vertex models ( 38 \ 
The specific character of the interaction leads to a definition of vertex states and related 
energies, which are different from those in the vertex models studied in the literature. The 
symmetries and the conservations typical of the solved classical model do not apply in our 
case, except for a special choice of the interaction constants. We show that, in the general 
case, the bare model can't have thermotropic phase transitions. 

In section 3. we take into account the conformational degrees of freedom and we show 
that the vertex states have an internal degeneracy. Therefore the model may exhibit a 
thermotropic transition of finite order. The thickness of the surface at critical point is also 
studied. 

In section 4. we analyze the critical dependence of the energy parameters of the model 
on the length of the molecule flexible section. We obtain the dependence of the order of 
the transition and of the critical temperature on such length. We then compare our results 
with some experimental data. 

In the last section we point out two lines of developement of our model: the first in 
view of a description of bilayers; the second toward an interpretation of the multiplicity of 



mesophases in liquid crystals. 



2. THE BARE MODEL 



We are interested in modelling a surface exhibiting orientationally ordered phases, so 
that we will assume the existence of a fixed director along which the molecules are aligned. 

At first, we don't consider the flexibility of molecular tails: in the bare model the 
conformational degrees of freedom are frozen, the molecules are rigid objects. 

Representing the surface through the centers of mass of the molecules, we consider 
their projections in the orthogonal plane to the director; we assume that such projections 
define a square Lx L lattice, formed by two intertwined sublattices A and B. Each lattice 
site in A (B) is occupied by a molecule with head dipole oriented up (down) with respect 
to the lattice plane; in this way we realize the head-to-tail setting of molecules, which, as 
described in the introduction, has been observed in smectic layers^ 11-14 ). 

The translational degrees of freedom of the molecules are constrained along the direc- 
tor; moreover, we don't consider vibrations (relative oscillations along the director, around 
the configurations which minimize the two-body energy). 

The surface can now be characterized by the scalar field of the altimetric heights of 

the molecular centers of mass, relative to a given reference. 

For a pair of molecules corresponding to adjacent sites in the lattice, the competition 



between dipole-dipole and chain-chain interactions determines two stable link states (with 
minimal energy). If the attractive chain-chain interaction prevails, we have two neighbour- 
ing molecules with tails facing one another and head dipoles forced in a non equilibrium 
(relative) position (see fig. la); this configuration defines the link state (0). If, instead, 
the antiferromagnetic interaction between head dipoles prevails, we have the link state (1), 
representing facing heads and tails forced to point toward opposite directions (see fig. la). 
Note that the direction of the line containing the centers of mass (local surface orientation) 
with respect to the director is different in the two configurations. 

Both link states (0) and (1) are minima for the pair interaction. In fact, a displacement 
around the configuration 1, increases the dipole-dipole contribution to interaction energy, 
whereas the tails are too far to compensate. Analogously, a displacement around 0, de- 
creases the attraction between the tails, but the opposite dipoles are again too far. The 
change from link-state (0) to (1) is obtained by shifting an A site molecule downward, or 
by shifting a B site molecule upward (fig. lb); since vibrations are excluded, the single 
molecule may only shift up or down of one unit length, relatively to its neighbours. 

In the bare model, given a distribution of links on the lattice, the height of a molecule 
in the site (I, J), relative to the height of the molecule in the site (0,0) (of type A), is a 
function of the values (0 and 1) of link states. We have: 

^j = E(-)^o 1 i + (-) J E(-)W a) 

j=0 i=0 

where Xij, yij are the values of the links + 1) and (i + 1, j) respectively. 

Note that the sign of the height variation between a site and its nearest neighbour depends 



on the parity of the site; furthermore, the sum of the height variations, over any loop on 
lattice sites, must be zero. This is a constraint for the distributions of links on the lattice. 

Let us define a plaquette as a loop of four links and a plaquette state as the set of 
values of such links, the constraint selects 6 of the 16 possible plaquette states, namely 
those where: 

xi,j-yi,j+i-x I+1 , J + y I ,j = (2) 

The previous description can be given in the formalism of vertex models ( S8 \ 
The plaquettes are identified with the sites of the dual lattice and the plaquette states with 
the vertices on such sites. A vertex is defined by four arrows entering or exiting the site 
(arrows-in/out) and connecting it with the nearest neighbours; it is easy to set a one-to-one 
correspondence between plaquette states and vertices. 

The so-called 8 vertex (8V) model is built on vertices with an even number of arrows-in 
(8 allowed states per site); the 6 vertex (6V) models have 6 allowed states per site. In 
particular the vertices of the classical 6V model must satisfy the so-called ice rule: two 
arrows-in and two arrows-out. 

The configurations of the model are sets of compatible nearest neighbouring vertices; the 
energy of a configuration is the sum of the energies assigned to the vertices. In the classical, 
exactly soluble models, the energy is invariant with arrows reversing; if such symmetry is 
broken, one has a model with external field. 

The allowed plaquette states of our model are mapped in 6 of the 8 vertices of the 
8V model, as represented in fig. 2.; this procedure formally leads to a 6V model, but, in 
comparison to the classical 6V, we have a different choice among the 8V vertices (vertices 



7 and 8 instead of 3 and 4). 

The configuration energy in the bare model is the sum of plaquette energies; for symmetry 
reasons the plaquette states 5,6,7,8 must have the same energy (e), whereas the energy of 
states 1 and 2 can be set to the values — uj/2 and +cu/2 respectively, by an appropriate 
choice of the zero for the energy scale. 

In the special case u = (vertices 1 and 2 with the same energy) , we can use the symmetry 
properties of the partition function of the 8V model with zero field (ref 38.ch.10) to prove 
that the partition function of the bare model is equal to the partition function of the classical 
6V model with ice rule^ 45 ^ 

Therefore, in this special case, our bare model is equivalent to the 6V model and it can be 
exactly solved. 

But, there is no physical reason to assume that tail-tail interactions have the same energy 
as dipole-dipole interactions, so that, in general, vertices 1 and 2 have unequal energies. 
In this case, the symmetry property, which is crucial in proving the above statement, 
cannot be generalized; it follows that, in a physically meaningful regime, our bare model 
is not equivalent to the classical 6V. Indeed, we have a specific version of the 8V model 
with external field (u); even if, to our knowledge^ 33 ' 38 ), the 8V model has not been solved 
exactly with non vanishing field, the specific features of the present model make it easy to 
study. 

From the convexity properties of the free energy and of the order parameter as functions 
of the field u, we can prove that in the bare model the free energy is differentiable with 
respect to u except, at most, in u = 0.( 45 ) 



It follows that we have no chance for thermotropic transitions in the fully oriented physical 
regime, as far as rigid molecules are considered. 

In the next section we shall show that the special case to = has a physical meaning 
if the internal degrees of freedom of the flexible molecules are considered. 
For this reason, we recall here the results of the classical 6V model that are relevant to 
our description. 

- According to the exact solution^ 38 ), one has two possibilities. First: if the energy e is 
strictly negative, the system is disordered; at any temperature the free energy is analytic 
and the path of the point representing the system in phase space lies in the disordered 
phase region. Second: if e > 0, for low temperatures the system is in the antiferroelectric 
phase and has a phase transition, of infinite order, at the critical temperature (3c = 

to the disordered phase. 

- The classical definition of order parameter, becomes, in our case: 

d 1 
m = -— lnZ| w=0 = (jp[ni - n 2 \) (3) 

namely the expectation of the difference between the numbers of plaquettes of type 1 (ni) 
and 2 (712). 

- The correlation length £, which is obtained from the behaviour of g^f^jj-, - for large 
distance of the sites (I, J) and (/', J'), is finite for (3 > (3c and diverges for (3 < (3c (ref 
38.ch.8). 

- The susceptibility \ associated with m, defined as \ = In Z , turns out to be finite 
for (3 > (3c and diverges for (3 < (3c ( 35 ' 39 ). 

In order to refer our bare model to known results, the altimetric description suggests 



the comparison with the BCSOS model ( 34 > 33 ). Therefore, we represent each plaquette by 
the mean height of the four constituent molecules and we identify the order parameter m 
in eq. (3) with the roughness of the surface. 

We carried on a numerical calculation, by transfer matrix method, for the bare model with 
u> = 0; in fact, in terms of our link variables the transfer matrix has an iterative structure 
which is very simple and easily computable ( 45 ). In particular, we checked the critical 
behaviour of the order parameter m and of the associated susceptibility x- 

Our results about the order parameter, reported in fig. 3, approach, at the critical 
temperature Tc, the infinite system asymptotic form, predicted by the exact solution 
(m ~ (T c - T)' 1 / 2 exp{-A(T c - T)" 1 /2}). The associated susceptibility (fig.4a) and its 
critical divergence with lattice size (xl(Tc) ~ L, see fig. 4b) show a good agreement with 
Montecarlo calculations on BCSOS (4 °). 

With reference to fig. 3, we note that the order parameter, as a function of temperature, 
is antisymmetric with u — > 0± (we shall indicate these functions with m + ((3) = —m~ (/?)). 
This fact, which comes from the invariance of the partition function of the bare model 
with respect to the exchange u <-> — u, will be crucial in the following. 

We are also interested in evaluating the dispersion 5 2 (h) of the height, related to the 
thickness of the surface^ 40 ' 41 ). According to the BCSOS model, S 2 (h) is expected to behave 
as ln£ for f3 > flc and to diverge as InL with lattice size L for f3 < /3c ( 35 ' 41 ). From a 
phenomenological point of view, we shall associate such a divergence with the disgregation 
of the surface (see sec. 3) 



3. SWITCHING ON THE INTERNAL DEGREES OF FREEDOM 



In this section, we shall consider the conformational degrees of freedom and prove 
that their switching on gives rise to a thermotropic transition of finite order. 

According to Flory's model of free chains, the single molecular tail is described in 
terms of conformational states. As we did for the case of alkanes lamellae^ 46 ) , we use a 
collective conformational variable A and write the chain partition function as a hierarchical 
sum: 

Z Tot = 1 + Z A=1 + Z A = 2 + ■■■ (4) 

Each term in the sum corresponds to a class of states characterized by conformations 
with a fixed transversal width A (with respect to the director). The first term (1) is the 
contribution of the all trans chains; the second term (with A = 1) comes from 1-kink or 
1-antikink conformations and from all conformations with sequences of alternating kink 
and antikink ; the A = 2 term collects all the contributions from conformations with twice 
the unitary width, and so on. 

To construct the partition function of the system, we have to define the structure of 
the interaction between molecules which have different tail conformations and belong to 
different plaquette configurations. 

We assume that the multiplicity of available conformational states depends, for each 
molecule, on the type of plaquette the molecule belongs to; the intermolecular interac- 
tion energies, instead, are assumed to be independent of the particular conformational 
state (provided it is available). This is argumented as follows. Let us consider a dressed 



plaquette of type 1 (see fig. 2; the molecules have facing tails and may form conformational 
defects); since the tails are strictly packed, they can only contain defects with small width 
A and the interaction is independent of the defect site in the chains. Considering now a 
dressed plaquette of type 2 (facing heads), the molecular tails are rather free and then may 
contain conformational defects with large A. 

From a general point of view, these assumptions represent a packing-flexibility constraint 
whose specific realization (see sec. 4) will not affect the general results discussed in this 
section. 

We write the partition function of the system in the form: 



where the first sum is over the plaquette states {Pi} of allowed configurations and the 
second sum is over the conformational states of the molecule in the lattice sites V, con- 
strained by the packing- flexibility condition with respect to the four adjacent plaquettes. 
The quantity h(Sy) is the Flory energy of the state Sy and h(P) is the plaquette energy. 
We can factorize the packing-flexibility constraint in such a way that we can sum over the 
conformational states and, using the geometry of the lattice, share the result among the 
plaquettes. The total contribution to the Z function from a single dressed plaquette takes 
the form: 




(5) 



{Pi} 



i s v\p} 



8 



e-W^Z^p + Z [2] d 2 , P + Z [0] $x,p} 



(6) 



X=5 



where Z^^Zp] are the partition functions of a single chain with states compatible with the 



plaquette of type 1 and 2 respectively and Z[ ] refers to the other configurations. 



In terms of the plaquette energies ±^ (states 1 and 2 ) and e (states 5,6,7,8) we have: 

Z(Ae,a/)= ^e^Llx-i^iufc. 

.e^EitT-? 111 ^]^^. 

One easily sees that this partition function has the same form as in the bare model: due to 
the conformational degrees of freedom, the bare interaction constants become the effective, 
temperature dependent, energy parameters: 



u = u + \\n^ (8a) 

P ^[2] 



E = e-I M |M ( ^ )1/2] (u) 

P Z {2] A [ 2 ] 

The quantity is a function of both the chain length N and with fixed N, it can be 

[2] 

written as a ratio of polynomials ( 46 - ) ; for f3 — > oo such ratio tends to 1 and for /3 — > it 
tends to the ratio between the number of constrained states (small A) and the number of 
free states (large A). 

The physical meaning of the additional terms in expressions (8a,b) - for brief f^' and 
/|q - is the difference between the free energy of a chain without and with steric con- 
straints. 

The function f^\/3) is negative, monotone, continuously differentiable with respect to (3; 
f[0\\P) nas analogous properties. 



From what precedes it follows that, with fixed N and provided Z^j 7^ Z[ 2 ], one inverse 
temperature (3* may exist, such that: 




= — u 



(9) 



For any /3 7^ /?* , since the energy ZJ is not equal to zero, with the analysis in section 2., 
we exclude the occurrence of any phase transition. 

If (3 = (3* , the plaquettes of type 1 and 2 have zero effective energy, whereas the others 
have positive effective energy e{(3*). It follows that, at (3 = (3*, the equivalence condition 
with the classical 6V model is verified. 

In order to study the behaviour of the free energy of the system e, uj) as f3 — > f3±, we 
consider the derivative of f((3,e,uj) with respect to (3. We have: 



As (3 — > the partial derivatives g= and -g= approach the values of the order parameters, 
calculated at /?*, of the bare model with vanishing external field (a; — > 0±); in particular, 
due to the antisymmetry of m, we can write: 



Now we compare the solution (3* of eq.(9) with the critical inverse temperature (3c of the 
classical 6V model. 

If (3* > (3c, we have: 




(10) 




m±(/r) 



(11) 



m+(/T) 



(12) 



namely, the first derivative of the free energy is discontinuous. 



In this case (3* is the critical inverse temperature of a first order thermotropic transition. 
The latent heat at the transition AH is the difference between the right and the left 
derivatives calculated at (3* in eq. (10); we have easily: 

df (N) 

AH = 2m((3*)[P*^-y=(3*] (13) 

g f (N) 

In the above expression, the quantity f3* ^ represents the melting heat of the tails; in 
fact, /jjj is the difference between the free energy of tails with straight conformations and 
with completely disordered conformations. Since the order of magnitude of m is 1, the 
transitional latent heat and the melting heat of the tails are really comparable. 
It is interesting to consider the critical behaviour of the following variance, which, in the 
limit of infinite size, represents the thickness of the surface: 

5 2 (s)=5 2 (h)+5 2 (D) (14) 

In this expression, the quantity S 2 (h) is the dispersion of the height of the molecular 
centers of mass, as defined in the bare model; the quantity S 2 (D) is the dispersion of 
the end-to-end distance of the molecular tails. At /?*, the dispersion S 2 (h) is finite and 
continuous (in fact, from the classical results, since (3* > (3c, one has S 2 (h) ~ ln£ with 
£ finite). On the other hand, since the molecular tails undergo conformational transition 
at (3* (associated with the change from plaquette of type 1 to plaquette of type 2), the 
quantity S 2 (D) abruptly changes from the characteristic value of straight conformations to 
that of disordered chains: both such values are finite. This proves that, at the first order 
transition, the thickness of the surface is discontinuous, but finite. 

Last, we recall that the orientation of the longitudinal axis of the molecules with respect 



to the surface plane changes in going from link (0) to link (1); as a consequence the mean 
orientation of the layered molecules is discontinuous at the transition. 
In conclusion, from a phenomenological point of view, the first order thermotropic transi- 
tion changes the (mean) geometry and the mutual positions of the molecules, but preserves 
the surface. This behaviour has been observed in smectic-smectic transitions; as said in 
the introduction, experiments show that molecules maintain the layered arrangement in 
the smectic phases, even if with different layer thickness and molecular overlapping (see for 
instance the Sai — Sa2 transition described in ref.13); it is also shown that the orientation 
of molecular axis with respect to the layer is different in the various smectic phases ( 9 > 10 ). 

We now consider the case (3* < /3c- 
One has: 

m + (J3*) = m-(J3*) = (15) 

which implies that the first derivative of the free energy is continuous at /?*; in this case 
no first order transition occurs. 

Clearly, the bare model critical behaviours are again maintained around (3* ; in particular, 
the susceptibility x diverges. If we consider the specific heat c v , at (3* we expect the same 
divergence as x ( m fact, c v ~ f^x)- This behaviour suggests the existence of a second 
order transition. From what said in sec. 2, the fact that the zero field condition uJ = is 
reached at (3* < /3c implies that the correlation length diverges at the critical point. In 
this situation, also the variance S 2 (s) diverges, due to the size divergence of 5 2 {h). 
We associate this transition with the disgregation of the surface, in fact the divergence 
of S 2 (s) means that the long range correlation between heights at different sites is lost. 



With reference to the phenomenology of nematics^ 9 ' 10 ), we note that our model behaves 
like a smectic layer of a system undergoing S-N transition. In fact, we are in a fully 
oriented regime and we go from a phase with strong long range correlation of heigths (the 
layered structure of smectics) to a phase with completely uncorrelated heigths (the nematic 
disorder in the orthogonal direction to the layer). 

We conclude this section noting that if Z^j = Z[ 2 ], namely if the small A class ex- 
tinguishes all the available chain conformations, we have no transitions at all (of course, 
since the orientational order is preserved in our analysis). 



4. CHAIN LENGTH DEPENDENCE 



In this section we study the solution of eq. (9) in function of chain length N. The 
obtained results depend on the qualitative properties of /, u, e as functions of iV; more pre- 
cisely, they are determined by the fact that such functions must have the same asymptotic 
behaviour for large N but have different scaling law for finite size effects. 

In order to study eq. (9) we first give the explicit choice for the packing-flexibility 
constraint and the associated representation of /[i] and /[qj. 

We assume that molecules in plaquette state 1 can substain only trans chain conformations; 
furthermore, we assume that molecules in plaquette state 2 may have any conformation. 
Consistently, we take: Z^j = 1 and Z[ 2 ] = Z^J t . As for Z[ ], the structure of plaquettes, 



with two link states (1), i.e. conformationally "free" chains, and two link states (0), i.e. 
packed chains, suggests Z[ ] = (Z^Z^])^ ■ 

(Alternatively, one may describe a system where at most A = 1 conformations are com- 
patible with plaquette of type 1; in this case: Z^\ = 1 + Z^}\) 

The dependence on A of /[i] and /[ ] is completely determined by Flory's theory. We 
consider the extremal cases: with A < 4, we have no conformational defects and then 
/[i] = /[o] = 0; with large A, the function takes the asymptotic form: 

f^\[3)=Nf([3)+ V (N,[3) (16) 

where / is independent of N and the function rj represents the finite size correction. It is 
rather easy to verify that such correction decreases exponentially with A as a consequence 
of the (mono)dimensionality of the chain; for example, one may extimate f^\(3) with 
transfer matrix method. Analogous considerations apply to /[ j. 

The problem is to model the dependence on A of the energy parameters. 
Contrary to the previous functions, such quantities are strongly affected by finite size 
effects and we are just interested in a range of lengths (A ~ 10 40) where such effects 
have to be considered. In order to justify our choices in modelling, we stress here some 
points. 

Let us consider the Lennard Jones interaction between two molecular chains. With faced 
and straight chains, as in the link state (0), the interaction energy behaves asymptotically 
like A x atan(N); this implies linearity in A, with power law corrections (— + • • •) for 
large A. With shifted and disordered chains, as in the link state (1), the interaction has, 
instead, a dominant term which is independent of A, again with power law corrections. 



Now consider the dipole-dipole interaction between molecular heads: the interaction energy 
is independent of N in the link state (1) (faced dipoles) and is high order in jy in the link 
state (0) (N measures the dipole-dipole distance). 

Given a plaquette configuration, the contribution of next nearest molecules in the plaquette 
is essential in determinig the sign of the total first correction in to the plaquette energy. 
For example, in the plaquette configuration 2, the next nearest neighbours chains are faced 
and statistically closer than in the plaquette configuration 1; this fact may change the sign 
of the jj term in the sum which determines u. 
Carring on such an analysis, we obtain the following general form: 



and the analogous for e. In this expression, Q(N) is the finite size correction, with power 
law dependence on jy. In particular, for what said above, the sign of the first correction 
is characteristic within a given system. 

We now search the solution of eq.9 in the form (3*(N) = (3^ + 5f3(N) ; we obtain two 
equations. The first one is independent of N: 



the solution fixes the order of critical temperatures within a given class of homologues; 



u = Nu (l + tl(N)) 



(17) 



^0 = /(Ax>) 



(18) 



in fact, the quantity = 



l 



is the observed value of transition temperature for large 



k B [3, 
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N. 



The second equation, regarding finite size corrections, is: 



6P(N) 



-u 



Q(N) 



(19) 



d/3 



'oo 



Since the prefactor is negative, critical temperatures may increase or decrease with N , 
depending on the sign ofQ(N). 



The comparison between (3c = n(Iq+ui ) ' an<1 leads to determine one critical 

length N*: we have C (N*) = (3*(N*). 

According to our results, for homologues with N > N* the first order thermotropic transi- 
tion occurs, whose main property is to preserve the layered structure as it happens in S — S 
transitions; for N < N* homologues, instead, the second order transition occurs, associated 
with surface disgregation as it happens for a smectic layer at the S — N transition. 

These general results agree with observations of phase transitions in liquid crystals; 
as said in the introduction, many data are available about homologous series of one-chain 
molecules, with typical regularities (see for insatnce ref.10). In particular, a phenomeno- 
logical expression, similar to the one used for phospholipids and bilayers ( 47 > 48 ) ? describes 
the observed S-S transition temperatures as a function of the number of Carbons per chain. 
We stress that both the characteristic polynomial dependence on chain length and the in- 
finite chain transition temperature, which appear in such expressions, in our case come 
directly from eqs. 18 and 19 which give the solution j3* of the model. 

As an example, we report, in Fig. 5, the fitting of critical temperatures of S-S transition 
in alkyl esters'- 49 -* and alogeno terminated compounds with alkyl chain ( 50 ); in the figure 
caption, we report the parameters of the fitting for each series of homologues. 



5. CONCLUSIONS 



The model presented in this work is a simple, exactly soluble model of surface which 
takes into account the flexibility of the constituent molecules. 

It proves the existence of transitions with different character and order, controlled by 
the length N of the flexible section of the molecule (a single linear chain). 
The character of the transition is determined by the dominant term in the expansion of 
the interaction energy; the dependence of critical temperatures on N is instead determined 
by the non-leading terms in the same expansion. 

The predicted transition leads, with increasing temperature, from a layered phase 
(smectic) to a different smectic or to a nematic phase, through conservation or, respectively, 
disgregation of the surface. 

In particular, if the surface is maintained, the transition is first order; it is associated with 
the conformational collapse of molecular tails and with a change in the mean orientation 
of molecules with respect to the surface plane. 

While at low temperature the surface is characterized by exposed molecular heads, which 
are polar and hydrophylic, in the high temperature phase the surface is formed by molecules 
with aligned heads and exposed, liquid-like chains, which are hydrophobic and weakly 
reactive. Regarding 3D interactions, the surface has a "soft" character; e.g., we expect a 
temperature dependent solubility: in the high temperature regime the solubility process 
has an higher energetic cost, associated with the loss of surface polarity. 

As a conclusion we point out two lines of developement of the present work. 



The first regards the reciprocal orientation of nearest neighbouring head dipoles. In 
the present model we explicitly require that adjacent head dipoles are opposite; this defines 
sites of type A and B in the lattice. If the requirement is removed, nearest molecules may 
get configurations with equally oriented head dipoles. The work in ref.47 points out the 
relevant role played by reversing molecular heads in liquid crystal transitions. For our 
molecules, with polar heads, the phases in which such reversed configurations are favoured 
must exhibit macroscopic polarity. We think that, at low temperature, the energetic cost 
of dipole reversing is not compensated by the entropic gain, whereas, with increasing 
temperature, reversed configurations may be favoured. 

From a technical point of view, to realize this improvement, the model should be changed 
by introducing a new link state and a new energy parameter. In such a situation the 
analytical treatment of the model will get much more involved. 

The second line regards a more accurate description of tails interaction. We think it 
would be interesting to consider the multiplicity of minima (decreasing in depth) which 
corresponds to relative shifts of adjacent chains of multiples of the intrachain periodicity 
(monomer- monomer distance). Clearly, the link state associated with each one of these 
minima represents two partially facing chains, having the active section of length / < N. 
If we consider only one of such further minima, we can repeat the present analysis noting 
that, in this case, we have oc / and to oc 1(1 +u(N)). In so doing we find that the latent 
heat is proportional to / and independent of N, while the behaviour of critical temperatures 
is again determined by jj expansion. These results agree with measured latent heats in 
S-S transitions of liquid crystals ( 52 ) which are rather constant, regardless of chain length. 



We could think that the interaction among distinct surfaces determines a mean separation 
and, as a consequence, favours a specific value for the active length /. 
We find technically difficult to achieve the solution of a model with all simultaneously active 
minima; on the other hand, the experimental evidence of a variety of rehentrant phases and 
the analysis of such phenomena carried on in ref.51, suggest the complexity of the problem. 
In this regard, we note that if, in our model, we considered the oscillations around the 
minima, we would not obtain a real modification of the present description; with harmonic 
oscillations, e.g., we could exactly sum over the associated degrees of freedom and obtain 
simply a redefinition of the energy parameters. 



Figure captions 

Fig. la 

The link state (0) (left) is a local minimum for the attractive interaction between 
molecular tails. The link state (1) (right) is a local minimum for the dipole-dipole inter- 
action between molecular heads. 

Fig.lb 

The exchange from link state (0) to (1) is obtained by shifting an A site molecule 
downward or by shifting a B site molecule upward. 
Fig.2 

Allowed plaquette states and associated link states: (0) = full line ; (1) = dotted 
line. Correspondence with vertices in the 8V model and related numbering, according to 
ref.(38). 

Fig.3 

Numerical calculation by transfer matrix of the order parameter as a function of 
temperature for different lattice sizes. Dotted line represents the exact behaviour (see 
text); in our calculation A = 7r 2 Tc/16y / 2. 

Fig.4 

a - Numerical calculation of susceptibility as a function of temperature from transfer 
matrix results for different lattice sizes. 

b - Size dependence of susceptibility at Tc in agreement with xl(Tc) ~ L (dotted 
line) . 



Fig.5 

Fitting of critical temperatures as a function of the number of Carbon atoms in the 
alkyl chain. The fitting has been obtained with: 0* = 1.47993 - + ^0493 _ 5^0055 

for observed^ 49 ) S-S transition temperatures of alkyl esters (•); f3* = 1.14883 — 0, °^ 027 and 
j3* = 1.12009 — °- 1] ^ 465 for observed^ 50 -* S-S transition of alogeno terminated aminoketone 
compounds (o -Br and LJ -J). 
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